Self-organized branching processes: Avalanche models with dissipation 
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We explore in the mean-field approximation the robustness with respect to dissipation of self- 
organized criticality in sandpile models. To this end, we generalize a recently introduced self- 
organized branching process, and show that the model self-organizes not into a critical state but 
rather into a subcritical state: when dissipation is present, the dynamical fixed point does not 
coincide with the critical point. Thus the level of dissipation acts as a relevant parameter in the 
renormalization-group sense. We study the model numerically and compute analytically the crit- 
ical exponents for the avalanche size and lifetime distributions and the scaling exponents for the 
corresponding cutoffs. 
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I. INTRODUCTION 



Many driven systems in nature respond to external 
perturbations by a hierarchy of avalanche events. This 
type of behavior is observed in magnetic systems , flux 
lines in superconductors fluid flow through porous 
media || , microfracturing processes j| , earthquakes || , 
and physiological phenomena [||. In these systems the 
distribution of avalanche amplitudes s decays as a power 
law, D(s) ~ s~ r , thus suggesting an analogy with criti- 
cal phenomena. Self-organized criticality (SOC) was pro- 
posed |Q as a possible framework to describe those phe- 
nomena. Power-law scaling would emerge spontaneously 
due to the dynamics, without the fine tuning of exter- 
nal parameters such as the temperature. Various models 
have been proposed with the aim of capturing the essen- 
tial features of avalanche dynamics and self-organization. 
In particular, sandpile models stimulated an intense ex- 
perimental HU, numerical JTo| , |TT| | and theoretical 12 h|] 
activity. 

As in the case of phase transitions, mean-field theory 
represents the simplest approach that gives a qualitative 
description of the system. Mean-field exponents for SOC 
models have been obtained in different ways [fi"5|-[2l| , but 
it turns out that their values (e.g., r = 3/2) are the same 
for all the models considered thus far. This fact can eas- 
ily be understood since the spreading of an avalanche in 
mean-field theory is a branching process |^2) because an 
avalanche can be described by a front of "non-interacting 
particles" that can either trigger subsequent activity or 
die out. The connection between branching processes 
and SOC has been investigated, and it has been proposed 
that the mean-field behavior of sandpile models can be 
described by a critical branching process |2^ |2^| . 

However, the nature of the self-organization was not 
addressed by the previous approaches. In fact the 
branching process is critical only for a given value of the 
branching probability, while in sandpile models there is 
no such tuning. Recently, we have introduced the "self- 
organized branching process" (SOBP) [^Tj, a mean-field 
model that allows one to clarify the mechanism of self- 
organization in sandpile models. Moreover, the SOBP 



model can be exactly mapped onto a two-state sandpile 
model in the limit d — > oo, where d is the dimension of 
the system. 

In experiments it can be difficult to determine whether 
the cutoff in the scaling is due to finite-size effects or due 
to the fact that the system is not at but rather only 
close to the critical point. In this respect, it is important 
to test the robustness of SOC behavior by understand- 
ing which perturbations destroy the critical properties of 
SOC model. 

It has been shown numerically p^-|30[| that the break- 
ing of the conservation of particle numbers leads to a 
characteristic size in the avalanche distributions. Here we 
generalize the SOBP in order to allow for dissipation and 
we show, in the mean-field approximation, how the sys- 
tem self-organizes in a sub-critical state. In other words, 
the degree of nonconservation is a relevant parameter in 
the renormalization group sense |yj . 

In section II we derive the SOBP from a dissipative 
sandpile model. In section III we study the approach to 
the critical state. The critical exponents are evaluated 
in section IV, and the results are verified numerically. 
Section V is devoted to conclusions. 



II. MODEL AND MEAN-FIELD THEORY 

Sandpile models arc cellular automata with an integer 
or continuous variable Zi (energy) associated with each 
site i of a d— dimensional lattice. At each time step the 
energy of a randomly chosen site is increased by some 
amount. When the energy on a site reaches a thresh- 
old z c the site becomes unstable and relaxes by transfer- 
ring its energy to its neighbors according to the specific 
rules of the model. In this way, a single relaxation can 
trigger other relaxations, leading to the formation of an 
avalanche. The boundary conditions are chosen to be 
open, so avalanches that reach the boundaries release en- 
ergy outside of the system. After a transient, the system 
reaches a steady state characterized by a balance between 
the input and the output of energy. 

Let us now consider a particular sandpile model: the 
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two-state model introduced by Manna [pl[ . Energy can 
take only two stable values Zi = (empty site) and Zi — 1 
(particle). When z% > 2 the site relaxes, z t — > Zi — 2, and 
the energy of two randomly chosen neighbors is increased 
by one. This rule conserves the energy, in this case the 
number of particles, during an avalanche and leads to a 
stationary critical state. 

Some degree of nonconservation can be introduced in 
the model by allowing for energy dissipation in a relax- 
ation event. In a continuous energy model this can be 
done by transferring to the neighboring sites only a frac- 
tion (1 — e) of the energy lost by the relaxing site p8| . 
In a discrete energy model, such as the Manna two-state 
model, one can introduce as the probability e that the 
two particles transferred by the relaxing site are annihi- 
lated For e = one recovers the original two-state 
model. 

Numerical simulations |2^,^9| show that the two ways 
of considering dissipation lead to the same effect: a char- 
acteristic length is introduced into the system and the 
criticality is lost. As a result, the avalanche size distri- 
bution decays not as a pure power law but rather as 



D(s)~ a - r h.(8/8 c ). 



(1) 



Here h s (x) is a cutoff function and the cutoff size scales 
as 

s c ~ <T V . (2) 

The size s is defined as the number of sites that relax 
in an avalanche. We define the avalanche lifetime T as 
the number of steps comprising an avalanche. The cor- 
responding distribution decays as 



D(T) ~ T- y h T {T/T c ), 



(3) 



where hx{x) is another cutoff function and T c is a cutoff 
that scales as 



(4) 



The cutoff or "scaling" functions h a (x) and hr(x) fall off 
exponentially for x ;§> 1. 

To construct the mean-field theory, we consider the 
model as d — > oo i.e. for an infinite dimensional lattice. 
When a particle is added to an arbitrary site, the site 
will relax if a particle was already present, which occurs 
with probability p — P(z — 1), the probability that the 
site is occupied. If a relaxation occurs, the two parti- 
cles are transferred with probability 1 — e to two of the 
infinitely many nearest neighbors, or they are dissipated 
with probability e. 

Since d — > oo implies that the lattice coordination 
number tends to infinity, the avalanche will never visit 
the same site twice, implying that each site that receives 
a particle from a neighbor relaxes with the same proba- 
bility. The avalanche process in the mean-field limit is a 
branching process. Moreover, we note that the branch- 
ing process can be described by the effective branching 
probability 



where p is the probability to create two new active sites. 
From the theory of branching processes p2| , we know 
that there is a critical value, p = 1/2, or 



1 



P = Pc 



2(1 -e)' 



(6) 



such that for p > p c the probability to have an infinite 
avalanche is non-zero, while for p < p c all avalanches are 
finite. The value p — p c corresponds to the critical case 
where avalanches are power law distributed. 

Boundary conditions are important for the process of 
self-organization. We can introduce the "boundary con- 
ditions" in the mean-field theory in a natural way by al- 
lowing for no more than n generations for each avalanche. 
We can view the evolution of a single avalanche of size 
s as taking place on a tree of N = 2 n+1 — 1 sites (see 
Fig. g). Note that we are not studying the model on a 



Bethe lattice 32 1; i.e., the branching structure we are 



discussing is not directly related to the geometry of the 
system. The number of generations n can, nevertheless, 
be thought of as some measure of the linear dimension of 
the system. If the avalanche reaches the boundary of the 
tree, we count the number of active sites a n (which in 
the sandpile language corresponds to the energy leaving 
the system) , and we expect that p decreases for the next 
avalanche. If, on the other hand, the avalanche stops be- 
fore reaching the boundary, then p will slightly increase. 



P(l"€) 




p = p(l- e) 



(5) 



FIG. 1. Schematic drawing of an avalanche in a system 
with a maximum of n — 3 avalanche generations correspond- 
ing to N = 2 n+1 — 1 = 15 sites. Each black site (•) can 
relax in three different ways: (i) with probability p(l — e) to 
two new black sites, (ii) with probability 1 — p the avalanche 
stops, and (iii) with probability pe two particles are dissipated 
at a black site, which then becomes a marked site (©), and 
the avalanche stops. The black sites are part of an avalanche 
of size s = 6, whereas the active sites at the boundary yield 
<Jz(p, t) = 2. The total number of "stopped" sites are p — 2, 
and there was one dissipation event such that k — 2. 
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To make the above statements quantitative, consider 
the evolution of the total number of particles M(f) in the 
system after each avalanche: 



M(t + 1) = M (t) + 1 - <r(p, t) - K (p, t). 



(7) 



Here a is the number of particles that leave the system 
from the boundaries and k is the number of particles lost 
by dissipation. Since M(t) — NP(z — 1) — Np, we 
obtain an evolution equation for the parameter p: 



p(t + l)=p(t) + 



1 - <r(p,t) - n(p, t) 
N 



(8) 



This equation reduces to the SOBP model [gTJ for the 
case of no dissipation [k = 0). The implications of Eq. (||) 
will be discussed in the following sections. 



III. SELF-ORGANIZATION: THE PROPERTIES 
OF THE STEADY STATE 

In order to characterize the steady state of the SOBP 
model, we rewrite Eq. (||) in terms of the average values 
of a and k indicated by angular brackets. The aver- 
age number of particles (a n ) leaving the system from the 
boundaries in a system of n generations is computed |2^j 
because of the recursive nature of the process: 



(a n (p,t)) = (2p(l-e)Y\ 



(9) 



The evaluation of the average number of particles dis- 
sipated during an avalanche is somewhat more involved. 
We can first relate the average value of k to the av- 
erage number of sites /i where an avalanche does not 
branch — either because of dissipation or because the site 
was empty (i.e., the avalanche stops), 



(«> = 2 </i> 



pe + 1 — p 



(10) 



The calculation of (k) then reduces to the calculation of 
(fi). If we denote by a m the number of active sites at 
generation to, then /i is given by 



n-l 
m=0 



(cTr, 



Pm+1 



l + s- 2a n 



(11) 



where s — J2m=o ° m ^ s the total size of the avalanche. 
The average value of s is obtained by summing the series: 



(s) = (vm) 



1 - (2p(l - e)) n+1 
1 - 2p(l - e) • 



(12) 



Combing Eqs. d)-©, one obtains that Eq. (^) in the 
continuum notation becomes 



dp 



N 



1 - (2p(l - e)Y 



pe + 1 - p 



1 - (2p(l - e)) 
1 - 2p(l - e) 



n+l 



2(2p(l - e)) n 



v(p,t) 

N 



(13) 



In Eq. (|l3|), we introduced the function rj(p,t) to de- 
scribe the fluctuations around the average values of a 
and k. We have shown numerically that the effect of this 
"noise" term is vanishingly small in the limit N — > oo. 

Without the noise term we can study the fixed points 
of Eq. (|l^). We find that there is only one fixed point, 

V = 1/2, (14) 

independent of the value of e; the corrections to this value 
are of the order 0(1 /N). By linearizing Eq. ([T3|), we find 
that the fixed point is attractive. This result implies that 
the SOBP model self-organizes into a state with p = p*. 
In Fig. H we show the value of p as a function of time 
for different values of the dissipation e. We find that in- 
dependent of the initial conditions after a transient p(t) 
reaches the self-organized steady-state described by the 
fixed point value p* = 1/2 and fluctuates around it with 
short-range correlations (of the order of one time unit). 
The fluctuations around the critical value decrease with 
the system size as 1/N. Thus it follows that in the limit 
N — > oo the distribution <j>(p) of p approaches a delta 
function 

0(p)~*(p-p*). (15) 

By comparing the fixed point value (|lj) with the criti- 
cal value (0), we obtain that in the presence of dissipation 



(e > 0) the self-organized steady-state of the system is 
subcritical. Fig. |3| is a schematic picture of the phase 
space of the model, including the line p — p c of critical 
behavior (^) and the line p — p* of fixed points (pl|). 
These two lines intersect only for e = 0. 
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FIG. 2. The value of the control parameter p(t) as a func- 
tion of time for a system with different levels of dissipation. 
After a transient, p(t) reaches its fixed-point value p* — 1/2 
and fluctuates around it with short-range time correlations. 
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B. Avalanche size distribution 



FIG. 3. Phase diagram for the SOBP model with dissipa- 
tion. The dashed line shows the fixed points p* — 1/2 of the 
dynamics, with the flow being indicated by the arrows. The 
solid line shows the critical points, cf. Eq. M. 



IV. CRITICAL EXPONENTS 

In this section, we study the critical properties of the 
model. In the limit n>l we obtain analytical results for 
the avalanche and lifetime distributions for any value of 
p, the effective branching probability defined in Eq. (|^). 
We show that the critical branching process with p = 1/2 
(obtained when e = 0) correctly reduces to the mean-field 
exponents r = 3/2 and y = 2. 



A. Generating functions 

The quantities P n (s,p) and Q n (cr,p) are defined to be 
the probabilities of an avalanche of size s and boundary 
size a respectively, in a system with n generations. The 
corresponding generating functions are defined by piJ 



f n (x,p) = ^2P n (s,p)x S 



(16a) 



by 



The solution of Eq. (|17a|) in the limit n 3> 1 is given 



f(x,p) 



1 - ^l-Ax 2 p{l-p) 
2xp 



(18) 



We expand Eq. ( jig ) as a series in x, and by comparing 
with the definition (16a), we obtain for sizes such that 
l«s<nH 



Pn( S ,P) = ^ 2{ \J )/7rP exp (s/s c (p)) . (19) 
The cutoff s c (p) is given by 



Sc(P) 



ln4p(l — p) 



(20) 



For avalanches with n < s < N it is possible to use a 
Tauberian theorem ]3^-|3f||, and show that P n (s,p) will 
decay exponentially. 

The next step is to calculate the avalanche distribution 
D(s) for the SOBP model. This can be calculated as the 
average value of P n (s,p) with respect to the probability 
density 4>{p), i.e., according to 



D(s) 



dp(j)(p)P n (s,p). 



(21) 



Since the simulation results show that cf>(p) for N ^> 1 
approaches the delta function S(p — p*) [cf. Eq. ([To])!, 
expression ( f2~l| ) reduces to 



D(s) = P n (s,p)\ p=pm(1 _ e) 

As a result we obtain the distribution 
^2"l + e-f 



D(s) = 



exp(-s/s c (e)) 



(22) 



(23) 



Jn(x,p) = ^2Q n {<J,p)3 



(16b) 



From the hierarchical structure of the branching process, 
it is possible to write down recursion relations for P n (s,p) 
and Q n (a, p) , from which we obtain p2| 



fn+l(x,p) = X [(1 -p) +pfn(x,p)] 



and 



9n+i{x,p) = (1 -p) +pgl(x,p), 
where fo(x,p) = go(x,p) = x. 



(17a) 



(17b) 



We can expand s c (p) in e with the result 



s c (e) ~ — , ¥ = 2. 



(24) 



Furthermore, the mean-field exponent for the critical 
branching process is obtained setting e = 0, i.e., 



3/2. 



(25) 



These results are in excellent agreement with the sim- 
ulation of D(s) for the SOBP model (cf. Fig. |). The 
deviations from the power-law behavior ( |23] ) are due to 
the fact that Eq. (Ill) is only valid for 1 < s < n Q. 
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FIG. 4. (a) Log-log plot of the avalanche distribution D(s) 
for different levels of dissipation. A line with slope r = 3/2 is 
plotted for reference, and it describes the behavior of the data 
for intermediate s values, cf. Eq. (jig]). For large s, the distri- 
butions fall off exponentially, (b) Data collapse produced by 
Eq. (Eh. 



FIG. 5. (a) Log-log plot of the lifetime distribution D(T) 
for different levels of dissipation. A line with slope y — 2 
is plotted for reference. Note the initial deviations from the 
power law for e — due to the strong corrections to scaling, 
cf. Eq. (||). (b) Data collapse produced by Eq. Q. 



C. Lifetime distribution 

The avalanche lifetime distribution D(T) is defined, 
for the model, as the probability to obtain an avalanche 
which spans m generations; here, we identify m with the 
time T. It follows that 



P(m,p) 



l+ i(cr = 0,p) 



(26) 



As for the avalanche distribution D(s) we have that 
D{T) = P(m — T,p) evaluated for p = p*(l - e). 
For p = 1/2 we use the general result p2] 



1 



l-Q ro ((T = 0,p) 
and obtain 



l + mp + OQnm), m > 1, (27) 



P(T,p) = (1 + OQxiT/T) + ...). (28) 

Note the strong correction to scaling to -D(T) in this case. 
Forp < 1/2 we find |§ 



l-Q m (fT = 0,p)~c 1 (2p) 1 



(29) 



for m 3> 1, where ci > is an unknown constant. This 
expression yields D(T) ~ eexp(— ne). 

We can combine the above results in the scaling form 



D(T)~T-«exp(-T/T c ), 



where 



V> = l. 



(30) 



(31) 



The lifetime exponent y was defined in Eq. (|^), where- 
from we obtain the mean-field result 



y = 2. 



(32) 



In Fig. we show lifetime distributions for different val- 
ues of e, together with the data collapse produced by 
Eq. ©. 



V. DISCUSSION AND CONCLUSIONS 

We have studied the effect of dissipation in the dynam- 
ics of the sandpile model in the mean-field limit (d — > oo). 
In this limit, the dynamics of an avalanche is described by 
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a branching process. We have derived an evolution equa- 
tion for the branching probability that generalizes the 
self-organized branching process (SOBP) introduced in 
Rcf. p7| . By analyzing this evolution equation, we have 
shown that there is a single attractive fixed point which 
in the presence of dissipation is not a critical point. The 
level of dissipation e therefore acts as a relevant param- 
eter for the SOBP model. We have determined analyti- 
cally the critical exponents describing the scaling of the 
characteristic size with e and the form of the avalanche 
distributions, and numerically verified the above results. 

These results prove, in the mean- field limit, that crit- 
icality in the sandpile model is lost when dissipation is 
present. It would be interesting to use a similar approach 
for other forms of perturbations. In particular it has been 
shown for other SOC models that the presence of a non- 
zero temperature p7fl or of a non-zero driving rate ]3§| ] 
are relevant perturbations leading to a non-critical steady 
state. 

Finally, we discuss the relations between the SOBP 
model and the simplest possible SOC system recently in- 
troduced by Flyvbjerg p9| . The minimal definition of 
SOC, as a medium in which externally driven distur- 
bances propagate leading to a stationary critical state, 
is well exemplified by the SOBP model. The disturbance 
is described by the branching process and the medium by 
the evolution equation for the density of particles in the 
system [Eq. (||)]. The example given by Flyvbjerg, being 
a two-state random-neighbor sandpile model, differs from 
the SOBP (27J in the way open boundary conditions are 
imposed. 
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